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We discuss the time spent by an electron propagating through a finite periodic system such as 
a semiconductor superlattice. The relation between dwell-time and phase-time is outlined. The 
envelopes of phase-time at maximum and minimum transmission are derived, and it is shown that 
the peaks and valleys of phase-time can be well described by parameters fitted at the extrema. For a 
many-period system this covers most of the allowed band. Comparison is made to direct numerical 
solutions of the time-dependent Schrodinger equation by Veenstra et al. {cond-mat/0411118} who 
compared systems with and without addition of an anti-reflection coating (ARC). With an ARC, 
the time delay is consistent with propagation at the Bloch velocity of the periodic system, which 
significantly reduces the time delay, in addition to increasing the transmissivity. 
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I. INTRODUCTION 



The question of "how long does it take for a non- 
relativistic particle to cross a barrier" , or in general 
anypotential, is a contentious one in quantum physics 
1, [j, El- There is for example the Hartmann effect 
J], according to which a transmitted particle might be 
found before it has reached the potential. Since 1990 
the subject has taken on renewed interest in the con- 
text of electrons in semiconductor superlattices (SL) 
I, i, 0, 1 S M EL HI 0, Hi HE IH. It is com- 
monly assumed that such propagation is ballistic, with a 
mean free path larger than the device dimensions. But 
because electrons can scatter from lattice vibrations, it 
is of interest to know how long the electron is exposed 
to such interactions. Many theories have been put for- 
ward as to how the time of passage should be defined 
and measured; among them dwell-time and phase-time 
are the best recognized. Veenstra et al. studied time 
dependence of propagation by direct numerical solution 
of the time-dependent Schrodinger equation (TDSE), us- 
ing gaussian incident wave packets. Their results agreed 
well with phase-time. In this contribution, based on a 
talk given at Theory-Canada 3 in June 2007, we explain 
why phase-time should be applicable in the context of 
semiconductor superlattices. 



II. TRANSFER MATRIX METHOD 

Assuming ballistic transport, a conduction band elec- 
tron in a potential cell of arbitrary shape on a < x < 
b satisfies the Schrodinger equation with a material- 



dependent effective mass: 



m*(x) 



dx 



dip(x) 
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ip{x) 



(1) 

The approximations leading to eqQ] are fully explained 
in Bastard's monograph [17]. At a layer boundary, il){x) 
and ip' = {h/m*)dtp/dx are continuous, to conserve flux. 
It is convenient to match ip( x ) 1° plane waves, normalized 
to unit flux, at the boundaries of the unit cell a < x < b 
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At either end, (B = L, R) is the wave number outside 
the periodic system, and Vb — tikis /nig is the velocity. 
We will consider systems without bias, so kr, — kjt — k, 
etc., but it will be convenient to retain the indices in order 
to know where various terms come from. The transfer 
matrix relates the coefficients on either side: 
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where 



M = 



l/t r* jt* 
r/t l/t* 



(3) 



in terms of the reflection and transmission amplitudes 
r(fc), t(k). It contains all the information about scatter- 
ing from that cell. One can show [20( that 



Ma 7 M* 



constant 



(4) 
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which says that the probability flux is preserved by the 
action of M, and which implies detM = 1. 

It is convenient to use Kard's parameterization [2l[ of 
M. In an allowed band we write 

Mn = cos0 — ism<fi cosh /i = 1/i = M 22 

M 2 i = -ie ix sin</> sinh fj, = r/t= M* 2 . (5) 

At a given energy, the Bloch phase is determined by 
TrM = 2 cos <fi; the impedance parameter from the ratio 
of |M2i|/ImMii = tanh/i, and the asymmetry parame- 
ter x from the phase of M21/M12 — exp (2«x). (x = for 
a reflection-symmetric potential cells. For the most part 
this paper is restricted to the symmetric case.) For a sim- 
ple square barrier-cell, is the ratio of average velocity 
outside to inside. Across an allowed band, <fi increases by 
7r, while \x is quite constant across most of the band, di- 
verging at the band edges. An example of this behaviour 
is shown in Fig. Q] taken from [l8j]. Panel (a) show cos0 
for a square barrier cell, and for a phase-shift equivalent 
gaussian barrier cell. Panel (b) shows the corresponding 
values for /1. The dotted lines in each panel are the pa- 
rameters of the single-layer ARC cell. The slope of cos 4>A 
is about half that of cos 0, while their Bragg points coin- 
cide. The rule of thumb for ARC's in optics would make 
HA = 



a) b) 




composed of cells, each described by M. Mathemati- 
cally, M is a hyperbolic rotation of a two-dimensional 
Dirac spinor, about an invariant axis whose polar angles 
are (,u, x) 0]. 

The usefulness of the transfer matrix for a periodic 
system arises from the property 

AfW = M N (<t>,n,x)=M(N4,,fi, X ) => 

— = cos N(j> — i sin N<t> cosh /i . (8) 

Hence the transmission probability for iV-cells 

|tjv| 2 = [1 + sin 2 N4> sinh V] -1 > V cosh 2 fi (9) 

shows narrow peaks determined by N<fi(E) = mir, m = 
1,2, ••■TV — 1. Between these peaks the |tjv| 2 touches 
the envelope of minima, l/cosh 2 /i. At the center of an 
allowed band, the envelope of minima crosses the curve 

In forbidden bands, the Bloch phase <j> — > pit + i6 
acquires an imaginary part; as a result the transmis- 
sion goes rapidly to zero. Pacher et al. [13] used this 
property to design an electron band-pass filter. Further, 
by adding a quarter-wave cell (anti-reflection coating, or 
ARC) at each end of the periodic array [H, [24]], they 
were able to increase the average transmission within the 
band from about 25% to about 75%. In a series of pa- 
pers 0, [H, [H, [13], some of us showed how to design 
an ARC which gives optimal transmission within a given 
miniband, by adding suitably configured potential cells 
on each end of a periodic array. Without an ARC, elec- 
trons which are transmitted do so via narrow resonances 
as in eq. [H so one expects a significant time delay. With 
an ARC, the incident plane wave is transformed into a 
Bloch wave of the periodic potential, and travels at the 
Bloch velocity. This should reduce the transit time. 



FIG. 1: Comparison of (a) Bloch phases of a square barrier 
cell (dash line) and a gaussian cell (solid line); also shown 
is cos <f>A of a single-cell ARC (dotted line) . (b) same for 
impedance parameters jj,, ha- 

In the Kard representation, transfer matrix factorizes: 

M = cos <f) — i sin <j> [cosh [i + sinh [i (cos x°x + sin x&y 

— g— <(x/2)o-» e (/i/2)cr x e -itj>a- z e ~(^/2)cr x e -i(x/2)a z 

The matrix of eigenstates is easily seen to be 

JJ = e -i(x/2)<T ie (M/2)cr x 

Substitution into eq. Ogives 



III. TIME DELAY 
A. Dwell time and phase time 

As already mentioned, the subject of time-delay in 
scattering is controversial. Razavy's book [3| is a useful 
a z introduction. Nussenzweig has argued that dwell- 

time is the best founded, for description of wave packet 
scattering. First we explain two of these theories. 

Classically, 



dt = 



dx 



dx 



v(x) hk(x)/m 
For a plane wave normalized to unit flux, as in eq. 



(10) 



MU = U e 



U 



;-«/ 2 cosh/i/2 e- lx / 2 sinh/i/2 
e l x/ 2 sinh^/2 e^/ 2 cosh///2 



1 , / „n,o m* 1 
\Mx,E)\ 2 = TT = — 



(11) 



(7) 



The eigenvalues are e ^ and e +l ^. Physically, the eigen- 
states are the Bloch waves of an infinite periodic array 



hk v c i 

Dwcll-time delay is defined to be the difference 

td = n [\i>{x,E)\ 2 - \Mx,E)\ 2 ] dx . (12) 
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Asymptotically, the scattering wave function, for waves 
incident from the left, is usually written 



problem (incident waves from left, or right). The S- 
matrix is symmetric and unitary, and may be written 



1 



\e tkLX +re~ lkLX ] 



1 



X < XL 
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\t\e ir > 
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i{2r/-5) 



(17) 



te 



ikax 



where 



t = \t\e"> 



x > x R 



(13) 



define the scatterin g ph ase shifts. 

Following Smith [29j], the Schrodinger equation gives 
the identity 



where f = —r*t/t* is the reflection amplitude for waves 
incident from the right. |r| 2 + |i| 2 = 1. 

Smith [29l | defined the time-delay matrix for a many- 
channel system as 



rl> = (H 



E) dE 



T = - iTl s^-- 

dE 

(\r\' + i\r\5')e lS 
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-ins j x 
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i\twy 



-flrf + i\r\(2r]' - 5'))e^-^ 



ip*ip 



2m 
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dx 



d_ 

dx 



<hp_ 
dE 



ip*ipdx 



h 2 

—Re 

2m 



r d_(d± 
V dx\dE 



dip dip* 
~dE~dx 

dip dip 



(18) 



dE dx 



where primes mean derivative with respect to energy. Af- 
ter some work we have 



XL 

(14) 



rii = +h |*| 



> dr] 
dE 



For a potential on a < x < b, xl < a and xr > b 
are positions at which a position measurement could be 
carried out. Taking only the scattering wave function in 
eq. [H 



T12 



T22 = 



dE 



dE 



1 dE J 



tan-' L ( J r f) +hr*t 



dr) 
dE 



1*1 
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dE 
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(19) 
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He 
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dE 
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H dE 
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In the case of a reflection symmetric potential, f = r, 
|^| which requires r] = 5 + it/2, so their derivatives are equal. 
As a result, all elements Ty are real, and the diagonal 
elements are equal. For a reflection-symmetric potential 
we have the compact result 



XL 
VL 



(1 



* dS \\ 



sm(2kLXL — S) 



(15) 



Til 



Tl2 



dr] 

n — = T22 
dE 

Ji-^r^ (sin -1 
dE V 



r\)=h(\r\'/\t\) 



(20) 



Using eqs. [TBI and [T4l in eq. [T2l gives the dwell-time as 



TD 



= h 



dE 11 + dE 1 1 



f5_f£)|t|» 



— — sm{2k L x L 



VR 



( 2x L 

\ V L 



8) 



(16) 



Since the potential lies between the limits, we can assume 
that Xl is a negative distance, while xr is positive. The 
terms on the second line have an obvious interpretation 
as the "free passage time" , so the top line is the "dwell- 
time delay" . But be careful: Nussenzweig @ has a lucid 
discussion of the case of an incident wave-packet, and 
finds an additional effect which arises from the uncer- 
tainty principle: you cannot localize a quantum particle 
in less than a de Broglie wave length. This effect cancels 
out between the free and interacting situations, so it af- 
fects dwell-time delay, but not dwell-time. We skip over 
this complication and rely on eq. [TBI 

Phase time was introduced by Wigner and Eisenbud 
[HI , and similarly relates time delay to the rate of change 
of the scattering phase. Scattering in ID is a two-channel 



The diagonal elements are called the "phase-time delay" . 
They agree with the first term in the top line of eq. 1161 
But due to coupling, the time associated with a process 
depends on just which mixture of the two channels is 
involved. The difference is that eq. [11] was derived as- 
suming that one has specified the channel with incident 
waves from the left, ignoring the other channel. The os- 
cillatory term of eq. [JB]has been subject of much debate 
also. Obviously it arises from interference between the 
incident and reflected waves. Winful [3l[ intepreted this 
term as the time taken to cross a distance of the or- 
der of the scattering length. On a more practical note, 
if xl — ► —00, and one averages fc^ over a wave packet 
narrow in energy, this oscillatory term is exponentially 
small. Therefore in a practical sense, either dwell-time 
or phase-time are equivalent for SL scattering, at least 
for symmetric potentials. 

Eq. [20] shows that there are two situations where the 
channel coupling vanishes: at a maximum or minimum of 
transmission. In a SL, the number of maxima is N—l, so 
as N increases, these points of uncoupling become closer 
together. 
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B. A tale of two phases 

In the transfer matrix method, the wave functions are 
defined with a different phase than in usual scattering 
theory: we reset the phase to zero on each side (x — a, b) 
of the potential array, rather than at an arbitrary origin. 
That makes M translation invariant, so M N describes a 
periodic system without having to include a shift in the 
origin. 

Adopting Nussenzweig's notation [ipL i iPr] for the 
asymptotic wave function to left and right [5j, we com- 
pare our ip(x) with the usual convention ip(x) as follows: 



ip(x) ~ [eW*-") +re -iH*-°) ■ te lk{x - b) ] 



\e lkx + fe %kx : 



te ikx 



e ika tp(x) ~ [e ikx + re -**(*-a») ; te ik{x ' w) ] , (21) 

where w — b — a is the total width of the potential. It 
follows that the phase r\ of our transmission amplitude 
is related to the usual phase by fj = i] — kw. Then the 
standard phase-time delay is 



Tph 



. drj 
l dE 



drj 

Je 



w 



dk 

Je 



drj 

Je 



w 

Vcl 



.(22) 



We conclude that the phase 77 of our transfer matrix am- 
plitudes gives phase time, not time-delay, because the 
free-passage time has to be subtracted from it. 



IV. PHASE TIME FOR SUPERLATTICE 
TRANSMISSION 

A. Relation of scattering phases to Kard 
parameters 

The S-matrix and the transfer matrix contain the scat- 
tering information in different forms. To see how phase- 
time applies to a SL, we examine the relation between 
the two descriptions. Equating 

— = cos Nip — i sin N<j) cosh fi = — l — e~ ir,N (23) 

we find the following relations between the two sets of 
parameters: 



cosr?Ar = |ijv| cos Nip 
tan t]n = tan Nip cosh fi 



sin^A? = |tjy| sin Nip cosh fi 

tan7?jv , tan 77 

= cosh fi = 



tan Nip 



tan i 



(24) 



where r\ = r}\ is the phase shift for a single cell. In the first 
allowed band, < ip < tt; tan 77 varies smoothly; cosh/i 
diverges at the band edge, like 1/ sin ip. Zeroes and poles 
of tan j]n coincide with those of tan Nip. These points are 
the scattering resonances N<p m — mir and the minima of 
transmission N<p p = (p+ 1/2)%. Near the poles, T]n lags 



behind Nip, then it has to catch up at the zeroes Ntp m . 
This is nicely illustrated in Fig. The steeper slope 
of /t]n near ip rn makes for a longer time delay at those 
energies. 



B. Phase time near the maxima of transmission 

In an infinite periodic array the electron would move 
at the Bloch velocity. Let the time to cross a cell of width 
d be tbi- Write tp/d = k, the pseudo-momentum. 



dE _ 1 dE 
dip d 8k 
dcP 
H dE 



dm* 

TBI ■ 



d tbi 



(25) 



In the following, the prime will mean d/dE. We can 
show using eqs. |2"31 and |2~41 that 

drjM [Nip' cosh/j + sin2iV</> sinh/x ///2] 

[1 + sinh 2 fi sin 2 Ncp] 

l + sin2Nip tanh/i (//' /{2N(p')] 



dE 

Tph = Nt B i cosh fi 



[l + sinh 2 fi sin 2 Nip] 



(26) 

Away from resonance, the denominator is a factor of 
|tjv| 2 which cuts off the phase time very sharply, causing 
it to mimic the shape of the transmission curve. Near a 
resonance, cp ~ mn/N + e, so that 

H m sin Nip ~ siniVe ~ N(E - E m )<p' m . 

This allows us to approximate | i jv | 2 as a Breit-Wigner 
resonance, with half-width r m = 2/(./V sinh /i m (p' m ): 



\tN\t 



1 



,E-E„ 



m \2 



r m /2 

Similarly, in the same vicinity, 
. dri N 



(27) 



dE 



l + 2b m ( 



N tbi. 771 cosh/i„ 



with 2 br, 



E — Err, 

r m /2 

1 



l + ( 



E-E, 



m \ 2 



N(p' m cosh/j r . 



r m /2 

1 

coth fj, m - 



(28) 

This is called a Fano resonance shape [3O] - From eq. 
[26l we see that the locus of phase time at transmission 
maxima is 



Tph.: 



Ntbi cosh/_i 



(29) 



Similarly, at transmission minima, sin 2Nif> — 1; the 
denominator of t p h (eq. [2^)1 becomes cosh 2 //, giving a 



downside locus of 



Ntbi I cosh /j 



(30) 



Tph^wTxw 

for phase time at transmission minima. The Bloch time 
for cells is the geometric mean of the two loci. 
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C. Phase time near the minima 

Transmission minima occur at cf> p = (p + 1/2)tt/N, for 
p = 1, 2, • • -n — 2. Close by, the denominator may be 
written 

tjj 2 = cosh 2 — tanh 2 fi cos 2 N<fi) 

~ cosh 2 fi p [l + ^ptanh/i p (<5i?p)] 2 [l — {N<fi p tanh^ p SE P Y 



cosh 



where 5E„ 



1 



M p ( SE p 



Nct>' p \Y p /2 



1 



SE p 
T p /2 



E - E p . The width 



r p = 2/[7V0 p tanh^ p ] (32) 
is large compared to the widths r m of the resonances. 



(31) 
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FIG. 2: Play model: coscfi, and angles 0, and rj, in units 
of 7r/2. Remarkably, cj> is reasonably linear over much of the 
band. 
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[ 1 + C P ( r p /2 )_ 




cosh n p 
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-2i?p(^)+(X?-l 





where C p = 



2 



and 



N<j>' 



(33) 



(34) 



Both C p and D p are of order 1/N. The prefactor is the 
locus of time-delay at minima: 



Ntbi 



NTi d(f> 



cosh fi cosh /j 



d cos d> , , 
= -Nh—^/lmM ll (35) 



which involves only well-behaved single-cell quantities. 

Veenstra et al. [18| compared their computed time de- 
lay to phase-time delay, and found good agreement. In 
particular they found that the locus of maxima and eqs. 
l27l l28l accounted for the overall picture. In an interest- 
ing paper, Pacher, Boxleitner and Gornik [3] pointed out 
that all theories of time-delay agree at the transmission 
maxima, so they concentrated their attention on those 
points and found the locus of maxima. The main differ- 
ence between their work and the present one is that they 
spoke of the mean velocity for an electron traversing the 
SL, rather than the time. Further they expressed their 
results in terms of the real and imaginary parts of the 
transfer matrix elements My, e.g. eq. I35[ rather than 
the Kard parameters. It is our opinion that the Kard 
parameters, given their simple behaviour, make the ex- 
pressions more easily understandable. 



D. Play Model 

To illustrate the above results, we take a simple model 
in which we specify cos <f> to be linear in energy in a band 
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FIG. 3: Play model: Transmission probability for one cell, 
nine cells, and envelope of transmission minima. 



between E = 50 and 75 meV. In addition, \t\ 2 is specified. 
All other parameters follow from these two. With such 
a model it is easy to see how changing some parameter 
will affect the results. 



cos<£ = 0.08(62.5 -E) 
\t\- 2 = 1 + 160/E . 



50 < E < 75 (meV), 



The phase 7/ is determined by cost? = 
impedance parameter follows from 

cosh/i = sin r\j (| t\ sin < 



(36) 



\t \ cos 6, and the 



(37) 



It diverges at each band edge, since |t| sin 77 is a smooth 
function of energy, while <fi runs from pir to (p + l)n. 

In Fig. [5] we show the play model cos = X(Eb — 
E), and the corresponding phases <fi and r\. r\ is quite 
linear, crossing <fi at the Bragg point. In this model, 
(j) 1 = A/sin0, with A = 0.08 meV^ 1 . 

In Fig. [3] we show the play model transmission for 
9 cells. The envelope of transmission minima is simply 
1/ cosh 2 fi, which gives direct physical significance to the 
impedance parameter. 



6 




E(meV) E (meV) 



FIG. 4: Play model: phase time, and loci of phase time at 
maxima and minima of transmission. Transmission probabil- 
ity is shown for orientation. 
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FIG. 5: Play model: t)n and N(j>, in units of n. Note that 
slope of t/n is largest at integer multiples of n. For clarity, 
only half the allowed band is shown. 



FIG. 6: Transmission probability compared to resonance ap- 
proximation (eq. I27[) . 
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FIG. 7: Play model: phase time, and its approximation by 
resonance formula eq. 1281 Locus of phase time at maximum 
transmission is shown to demonstrate that phase-time max- 
ima are very little shifted from the locus. 



In Fig. 2] we show the play model phase-time, along 
with the envelopes of maxima and minima. In the back- 
ground for orientation are the transmission peaks, which 
line up well with the maxima of phase-time. 

In Fig. [5] we show the dance of the Bloch phase for TV- 
cells, and the corresponding transmission phase ry^r, in 
the lower half of the allowed band. The curve in the up- 
per half of the band is a double reflection of this, ending 
at 9ir. Except at the band edges, the lines cross at every 
half-integer multiple of 7r. Since r/N is catching up at in- 
teger multiples of n, the steeper slope leads to a longer 
phase-time at the transmission resonances. 

In Fig. [5] we show (solid line) the transmission and 
(dashed line) the Breit-Wigner resonance fitted at the 
peaks. This is truncated at two standard deviations; the 
horizontal lines simply connect the B-W curves between 
successive peaks. The agreement is excellent. 

Fig. [7] is similar to Fig. O but for the phase-time. The 
Fano-shape formula fitted to the maxima also does an 
excellent job of reproducing the exact calculation. 



E. Realistic potential model 

Some results using a realistic semiconductor potential 
of Pacher and Gornik [23|, [gH are shown in Figs. [SJ [9] 
Here we took five potential cells, giving four resonances 
in the band, a number chosen for comparison with the 
calculations of Veenstra et al. [Hj]. Fig. [5] corresponds to 
Fig. Q]for the phase-time. For this potential the envelope 
of minima is lowermeaning a larger fi. Also shown is the 
Bloch time, the geometric mean of the two envelopes. 

Fig. [5] can be compared with Fig. [7J The Fano for- 
mula fitted at the maxima of transmission agrees with 
the exact result over two standard deviations. In Fig. [S] 
we have included eq. 1331 fitted at the minima of trans- 
mission. The fit to the minima is not so good, largely 
because the half-width at minima is so much larger than 
at maxima. Still, we can say that between them the ap- 
proximations reproduce the phase time over about 90% 
of the band. If there were more layers, the peaks and 
valleys would be narrower and the agreement would im- 
prove, as in Fig. [JJ For large N the phase time can be 
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FIG. 8: Pacher-Gornik 5-cell array: phase time (solid line), 
along with the loci at maxima and minima (long dashes), and 
Bloch time (dash - dot), which is their geometric mean. The 
chain-line is four times the transmission probability. 
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FIG. 10: Time-dependent numerical solution of wave equa- 
tion for a Pacher-Gornik 5-cell array plus ARC: phase-time 
compared to Bloch time. Transmission shown in background. 
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fit at the peaks is excellent over two half-widths, but not so 
good at the phase-time minima. The long dash and dash-dot 
lines are the loci of phase-time at maximum and minimum 
transmission. 



calculated reliably using only properties at the extrema, 
where the various theories of time-delay agree with each 
other. 

Finally, in Fig. [10] we show the results calculated by 
Veenstra et al. for the Pacher five cell array, plus a two 
layer ARC. The dotted line shows the transmission is 
close to 100% over most of the band width. There are two 
lines shown for the Bloch time. One takes into account 
only the periodic 5-cell system, and the other includes 
a correction for the ARC layers on each end. The two 
estimates are close together. The time delay extracted 
from the TDSE calculations is indeed very close to the 
Bloch time, which is what we expect from the argument 
that the ARC converts the incident plane wave into a 
Bloch state. 



V. CONCLUSION 

We have outlined the transfer matrix method for trans- 
mission of electrons in a one-dimensional AlGaAs/GaAs 
superlattice, using the Kard parameterization. We then 
compared the two most commonly used theories for the 
average time spent in crossing a potential region, namely 
phase time and dwell time. These were applied to SL 
transmission, first for a play model and then for the po- 
tential that corresponds to experiments of Pacher and 
Gornik. We noted that the phase-time is well defined 
both at maxima and minima of transmission. In the 
neighbourhood of these points, the transmission can be 
described by a Brcit-Wigner resonance formula, with the 
parameters extracted at the extrema. The same holds 
for the phase-time, except it is a Fano shape resonance. 
The fits at maxima are good over two half-widths, and at 
the minima over one. Taken together, this covers almost 
all the band width, for a system with more than a very 
few periods. 

As the number of periods N of the SL is increased, 
the widths of the peaks and valleys decrease, and the ap- 
proximate forms become more and more accurate. This 
explains the success of phase-time for describing the time 
spent in traversing a SL. 

The locus of phase-time at maxima was derived in 
Veenstra et al. [3, and by Pacher et al. [IJ]. The 
locus of phase-time at mimina is new. Their geometric 
average is the Bloch time for traversing a SL. An ARC 
works by converting the incident plane wave into a Bloch 
wave of the periodic system. Veenstra et al. [l8| studied 
the time delay for the SL plus ARC by direct numerical 
solution of the TDSE using gaussian wave packets. In 
Fig. [TU1 taken from that reference, it can be seen that 
indeed the time-delay with ARC agrees quite well with 
the Bloch time. The ARC not only increases the average 
transmission, but it smooths out the dwell time, remov- 
ing the peaks and valleys associated with the exponential 
decay of the resonances. 
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